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Синтез адаптивної системи керування 
електровибуховим перетворенням енергії 
на основі інверсної Шш7гу-моделі 


Побудовано інверсну модель електровибухового перетворення енергії на базі експериментальних 
даних з використанням методу їш77у-апроксимації, основаного на теорії нечітких множин. На базі 
інверсної моделі синтезовано нечіткі регулятори для корекції параметрів адаптивної системи керування 
електровибуховим перетворенням енергії в умовах непередбачуваних зовнішніх впливів. 


Вступ 


Результати теоретичного і експериментального дослідження електровибухового 
перетворення енергії в каналі високовольтного розряду в рідині як об'єкта керування 
показали, що об'єкт відноситься до дискретно-неперервних багатовимірних нестаціонар- 
них стохастичних систем |1). Оператор відповідності між координатами вхідного 1 
вихідного векторів об'єкта та статистичні характеристики інформаційних координат 
залежать від положення в просторі станів, тому для реалізації основної переваги 
електровибухового перетворення енергії - керованості, необхідно синтезувати адаптивне 
керування, яке забезпечує оптимальні режими в умовах нестаціонарності параметрів 
середовища та невизначеності зовнішніх збурень. 

Мета роботи - синтез нечітких регуляторів (НР) адаптивної системи керування 
на основі інверсної моделі електровибухового перетворення енергії, побудованої на 
базі експериментальних даних із застосуванням сучасних методів Шггу-апроксимації, 
що базуються на теорії нечітких множин. Результатом чого стане розширення зони 
керованості, підвищення точності підтримки оптимальних технологічних режимів в 
реальних умовах. 


Побудова інверсної Г/гу-моделі для синтезу НР 


Аналіз існуючих моделей та систем керування електровибуховим перетворенням 
енергії показав, що вони побудовані на основі лінеаризації об'єкта і забезпечують 
керування тільки в околі точки номінального режиму при певних припущеннях 
відносно збурень. 

Забезпечення керованості електровибухового перетворення енергії у всьому 
просторі станів потребує введення механізмів адаптації для врахування нестаціонарності 
та стохастичності об'єкта. Механізм адаптації можна реалізувати, використовуючи 
системи керування на основі інверсних моделей, які генерують реакцію на збурення, 
визначаючи керуючий вплив, що гарантує необхідні значення вихідних координат. 
Достовірність інверсної моделі забезпечується точністю та повнотою опису властивостей 
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об'єкта та його статистичних характеристик. Побудова інверсної моделі передбачає 
отримання зворотної функції об'єкта керування, тобто функції, що дає змогу знаходити 
таке значення керуючого впливу, яке забезпечить задане значення вихідної інформа- 
ційної координати, а відповідно, і заданий режим електровибухового перетворення 
енергії та необхідний технологічний вплив, при поточному стані контрольованих 
параметрів об'єкта. 

Аналітично пряму та інверсну модель електровибухового перетворення енергії, 
які адекватно описують об'єкт у всьому просторі станів, побудувати складно, так як 
фізичні процеси, що відбуваються в каналі розряду, неоднозначні, недостатньо 
вивчені, погано піддаються формалізації |2|, тому для отримання цих залежностей 
доцільно застосувати сучасні методи їш7гу-апроксимації на основі нечіткої логіки, 
використовуючи бази експериментальних даних. З метою отримання бази даних 
проведено теоретичне і експериментальне дослідження процесу дискретного, з 
періодом Т, електровибухового перетворення енергії в розрядному контурі генератора 
імпульсних струмів як об'єкта керування У(пТ) - АХ(птТ) (рис. 1). 


Хипт) У(птТ, т) 
цу А із (9- 
С(птТ) 


иг(б алі (7). 


циті РОЇ род 14-14 
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Рисунок 1 - Структурна схема електровибухового перетворення енергії 
як об'єкта керування 


Координатами вхідного вектора Х(пТ) можна вважати параметри розрядного 
контуру: У - зарядна напруга накопичувача; С - ємність батареї конденсаторів; 
І - індуктивність розрядного контуру, /(пТ) - величина розрядного проміжку та 
ріп| - питомий опір рідини, п - поточна реалізація процесу. Вихідний вектор 
У(пТ,т) - результат електровибухового перетворення енергії в каналі розряду, харак- 
теризується випадковими імпульсними функціями: (т) - розрядного струму, и(т) - 
напруги на розрядному проміжку, р(т) - тиску в каналі розряду, т - тривалість розряду 
(те «Т), Як координати вихідного вектора У(пТ) можуть бути прийняті функціонали 
цих функцій, які кількісно і однозначно їх характеризують, наприклад, максимальні 
значення розрядного струму іліп| 1 тиску ри|п| та пробивна напруга шпріп|. Експери- 
ментальне дослідження координат вихідного вектора |2| показало, що вони є дискретними 
випадковими функціями з нормальним законом розподілу в кожній точці факторного 
простору (відповідно, критерію згоди У, це не суперечить істині, з рівнем значимості 
а. - 0,01) ії можуть бути представлені у вигляді: у,|п)| - мулу, де М/|п| - математичне 
сподівання вихідної координати, у» -- завада - стаціонарна випадкова величина з законом 
розподілу Гауса, статистичні характеристики якої визначаються процесами формування 
каналу розряду і залежать від положення об'єкта в просторі станів. Оператор об'єкта 
А множини вхідних станів Х ставить у відповідність множину У функціонального 
простору можливих реалізацій вихідних функцій. Як статистична ефективна інформа- 
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ційна координата прийнята лінійна комбінація корельованих з від'ємним коефіцієн- 
том кореляції вихідних координат: Миі| з ішіп| Я Кипрі п|/йв| п|, що має мінімальну 
дисперсію. 

Інверсну модель будуємо у вигляді | п| - К(Мпі, д|лі|), для чого необхідна база 
навчальних даних, тобто множини груп сигналів (Мп), р|пі|), що будуть поступати на 
вхід системи керування, та відповідно очікувані (еталонні) значення керуючого 
сигналу | п|. Такі групи сигналів формуємо, використовуючи базу експериментальних 
даних, наведену в |3|. 

Першим кроком при побудові нечіткої інверсної моделі є розподіл простору 
вхідних і вихідних сигналів на діапазони, які будуть слугувати термами при побудові 
функцій приналежності. Використано прямий метод побудови функцій приналежності, 
який не вимагає абсолютно точного визначення функцій приналежності, а лише тип 
функції та характерні значення, в даному випадку - значення в експериментальних 
точках за планом експерименту, які є опорними для діапазонів, прийнятих в якості 
терми. Приймаємо трикутну форму функцій приналежності. Одна з вершин ФІ зна- 
ходиться в центрі кожного числового діапазону значень координати і їй відповідає 
значення функції и'(х), що дорівнює І, дві інші вершини лежать в центрах сусідніх 
діапазонів, їм відповідає значення функції, що дорівнює 0). Як лінгвістичні вхідні змінні 
приймаємо координати | пі|, 5 п). Кількість термів (лінгвістичних значень, а у нашому 
випадку термів - числових діапазонів) для кожної змінної вибираємо рівною кількості 
рівнів значень за планом факторного експерименту. В даному випадку кількість термів 
дорівнює 5 та 4. Функції приналежності (ФП -- т/; - ц'(х)) координат вектора стану 
ріп|, У |п| апроксимуємо трикутною функцією (рис. 2). Для вихідної координати, тобто 
керуючого сигналу || п|, виконуємо такі ж процедури: розподіл простору значень сигналу 
на діапазони, які слугують термами при побудові функцій приналежності; побудова три- 
кутних функцій приналежності, координати вершин яких є еталонними очікуваними 
значеннями вихідної координати | п| (рис. 2). 
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Рисунок 2 - Функції приналежності |(|п|, У; |п|, Пп) 


Процедура їиг7гу-фікації вхідних сигналів, тобто визначення приналежності по- 
точного значення координат вхідного вектора до того чи іншого діапазону, проводи- 
ться за максимальним значенням функції приналежності и'(х). 

База правил нечіткої інверсної моделі формується на основі бази навчальних 
даних, в якості якої використано експериментальні дані, за алгоритмом побудови бази 
нечітких правил на основі чисельних даних |4). До бази правил включається правило, 
яке має найвищий ступінь істинності. Таким чином розв'язується проблема супереч- 
ливих правил, а також значно зменшується їх загальна кількість. Блок-схема алгорит- 
му побудови бази правил на основі чисельних даних наведена на рис. 3. 
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Початок 


и 5В(Пу» 
5в(В,р) 
Так 
5 (К,ру- 
5 (Пп) 


Так 


Г Розділити простір значень Х, р, І на 
діапазони В, р, В. Побудувати 
відповідні їм функції приналежності. 


Створити таблицю БП для запису бази 
правил і таблицю 5В ступеню істинності 
| правил; заповнити ці таблиці нулями 


Ініціювати лічильник одиницею 


Вибрати наступну групу даних 
ХЙ, РІЙ, ПЛ 


Встановити ступені приналежності даних 
до діапазонів (нечітких множин) і 
сформулювати відповідне правило 

П: Якщо (2 В) і(рР), то (І-В)) 


" Встановити ступінь істинності 
правила П за формулою: | 
| 58 (П) з м ?/ (сіє м У (Уд' м У (І) 


( Вписати правило П в таблиці БП та 85В: 
5 (В,Рру- 5в(П) 


Інкриментувати лічильник 


Чи розглянуті усі групи даних 


Рисунок 3 - Блок-схема алгоритму побудови інверсної моделі 
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Після застосування наведеного алгоритму до бази даних, отриманих при експери- 
ментальному дослідженні електровибухового перетворення енергії (31, отримано несупе- 
речливу базу знань (табл. 1), яка є табличною формою запису бази правил нечіткого 
регулятора НРІ на основі нечіткої інверсної моделі. Значення інформаційної координати 
Мп| дано у відносних одиницях, за базове значення прийнято амплітудне значення розряд- 
ного струму при короткому замиканні Іа 


Таблиця 1 - Таблична форма бази правил НРІ 


Му р.Омм 

6,0 ЛЗ 10,0 15,0 20,0 
0,40 0,047 0,050 0,056 0,094 0,100 
0,45 0,040 0,042 0,048 0,081 0,090 
0,50 0,032 0,035 0,038 0,069 0,082 
0,95 0,025 0,027 0,030 0,052 0,065 
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База правил у вигляді нечітких логічних рівнянь дозволяє пов'язати вихідний сиг- 
нал Дп)| та ФП координат вектора стану (п), У|п|, в результаті чого отримуємо значення 
керуючого впливу для певних діапазонів числових значень координат вектора стану. 

Базу правил нечіткого регулятора НР2, що відтворює залежність середньоквад- 
ратичного відхилення інформаційної координати від положення в просторі станів сх - 
- Е (Цлі|, д|п)), синтезовано на основі експериментальної бази даних (табл. 2). Зна- 
чення середньоквадратичного відхилення сх приведене до математичного сподіван- 


ня 2, (лі). 


Таблиця 2 - Таблична форма бази правил НР2 


р, Ом м 

І, м 6,0 9 10,0 15,0 20,0 
0,025 0,028 0,024 0,023 0,015 0,010 
0,050 0,037 0,036 0,036 0,014 0,012 
0,075 0,042 0,040 0,039 0,015 0,013 
0, 100 0,074 0,069 0,044 0,017 0,016 
0,125 0,085 0,072 0,055 0,032 0,024 


Моделювання нечітких регуляторів НРІ, НР2 на основі отриманої бази правил 
проводимо в середовищі МАТІ АВ, використовуючи пакет Ки77у Ібіс Тооібох (51, 
який має простий інтерфейс для проектування 1 діагностики нечітких моделей. Графічні 
засоби Ки7гу Іобіс Тооїбох дають змогу інтерактивно відслідковувати поведінку 
системи. Результати моделювання: НРІ - /п| - Е(Мпі|, д|п|) та НР2 - сх- Е (пі), д|п1) 
за допомогою їш7гу-апроксиматора представлені на рис. 4, 5. 


Рисунок 4 - Модель НРІ: Ди|.Е(М пі, рГп|) Рисунок 5 - Модель НР2: с» -Е (Цпі|, дл |) 


параметрів в залежності від положення в просторі станів показана на рис. 6. Важливою 
властивістю системи керування з використанням НР на базі інверсної моделі є те, що вона 
забезпечує необхідні режими при будь-якому положенні об'єкта в просторі станів. Синте- 
зована система керування реалізована на базі мікроконтролера архітектури АУБ. Ресурси 
мікроконтролера дозволяють організувати досить складну ієрархічну багатомірну систе- 
му керування електровибуховим перетворенням енергії, що здатна вирішувати завдання 
керування в реальному масштабі часу із застосуванням різних алгоритмів керування у ви- 
гляді підпрограм. 
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5 - стохастичні впливи 
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Регулятор 


Рисунок 6 - Блок-схема адаптивної системи керування з НРІ, 2 


Висновки 


Використання інверсної моделі, побудованої на принципах нечіткої логіки, при 
синтезі нечітких регуляторів для корекції параметрів адаптивної системи керування 
електровибуховим перетворенням енергії в залежності від положення в просторі ста- 
нів забезпечить адаптивність системи керування при зміні технологічних параметрів та 
параметрів середовища в широкому діапазоні. Така система керування зможе генеру- 
вати адекватну реакцію на збурення за рахунок визначення відповідного значення керую- 
чого впливу за допомогою інверсної моделі. 

Синтезована система керування забезпечує ефективне керування електровибуховим 
перетворенням енергії в усьому просторі станів при невизначеності зовнішніх збурень 
та нестаціонарності параметрів середовища. 
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С.С. Козьрев 

Синтез адаптивной системьг управления злектровзрьвоопасньюм превращением знергии на основе 
инверсной їш/у-модели 

Построена инверсная модель злектровзрьтвного превращения знергии на базе зкспериментальньгх данньтх с 
использованием метода їшиг7у-аппроксимации, основанного на теории нечетких множеств. На базе инверсной 
модели синтезированьт нечеткие регуляторьт для коррекции параметров адаптивной системьт управления 
злектровзрьтвоопасньтюм превращением знергии в условиях непредсказуемьїх внешних влияний. 


5.5. Косупіеу 

ЗупіНпезіз ої АЧаріїуе Сопітої Зузіет їог Дізсрагое Кпегоу Сопуегзіоп, Ваз5еа оп Шуег5е Киггу Модеї 

ТЬе Гаггу шуег5е тодсі ої сопігої ої дїзспагее епегеу сопуєгзіоп ууа5 деуеЇоредй. ТПе плодсі 15 Базед оп 
ехрегітепіа! ага мів Га7гу арргохітайоп. Тре Ги77у іпуєг5е плоде! ргоуідез адаріабійсу ої сопігої 
зузієт шпаєг уагіабіе іесппоїоріса! рагатеїегя апа ехіегпа! сопантоп5. 


Стаття надійшла до редакції 14.05.2009. 
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